Role of local duality invariance in axion electrodynamics of topological insulators 
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Advances in material technology and confluence of ideas from particle physics, quantum field 
theory and condensed matter physics have led to the discovery of new states of matter as well as 
new physical phenomena: one of them termed as topological insulator has attracted a great deal of 
attention recently. Speculations on the possibility of observing the most elusive objects like axions 
and monopoles in topological insulators have led to studies that emphasize the role of symmetry 
and universality. In this paper we argue that electric-magnetic duality could be of deep significance 
in this context. We develop a duality invariant theory of topological insulators and show that under 
appropriate conditions this theory reduces to the axion electrodynamics for static case, topologi- 
cal quantization is related with the multi-valuedness of the duality gauge potential, and modifies 
Faraday's law for dynamical axion that would change the dispersion relation of axionic polariton. 
A new effect dual to the dynamical axion field effect is predicted and its physical consequences are 
discussed. 
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Spin-orbit coupling arises naturally in Dirac's relativistic equation of electron in external electromagnetic (EM) 
field [1]; recent advances in the study of new materials have shown its renewed significance in condensed matter. 
One of the striking examples is the spin Hall effect predicted theoretically and soon established experimentally [2, 3]. 
Spin Hall effect at a basic level is accumulation of the opposite spin states of electrons at the opposite edges of a 
material under suitable conditions of the application of EM field. Similar to quantum Hall effect it also gets quantized 
though the origin of quantization is different. Confluence of topology, symmetry (time reversal) and microphysical 
interaction mechanism (spin-orbit coupling) paved the path for the realization of so called 'exotic state of matter', 
namely topological insulator [2, 3]. A simple way to understand topological insulator is to consider its response to 
the applied EM fields, first recalling the standard Maxwell theory in which for a linear isotropic medium the response 
at a macroscopic level is embodied in the constitutive relations 

D = E + 4ttP = eE (1) 



H = B - 47rM = — (2) 
M 

Here P and M are electric and magnetic polarization vectors, and the scalar quantities e and /i are the usual 
permittivity and permeability respectively. A topological insulator (TI) is characterized by an additional parameter 
9 that couples to the pseudo-scalar product of electric and magnetic fields in the effective topological action [2] 
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and modifies the constitutive relations 



Stop = / —aE.B <?x (3) 
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D = eE -qB (4) 
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Here a = e 2 /Tic. A constant parameter 9 implies that the integrand in (3) becomes a total derivative, therefore Stop 
would not alter the field equations. 

Note that the term E.B in (3) transforms as a pseudo-scalar under space reflection (P) and it is odd under time 
reversal (T); since the standard action has to be a Lorentz scalar 9 must be odd under both P and T. Evidently 
in common insulators 9 — 0. Interesting possibility assuming T-invariance exists when 9 = 7T, and this defines 
the new class of insulators - TI. Authors in [2] emphasize that the parameter 9 corresponds to magneto-electric 
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polarization in their theory and expressions (4) and (5) embody the topological magneto-electric effect (TME) in 
which magnetic (electric) field could induce electric(magnetic) polarization. Another new result is that an electric 
charge near the topological interface induces a mirror or image magnetic monopole [4]. Intriguing generalization [5] 
is to break T symmetry via magnetic fluctuations so that 9 becomes a dynamical field variable - these are called 
topological magnetic insulators (TMI). 

Promoting 9 to a dynamical field entails that the integrand in action (3) would change the field equation since it 
cannot be written in the form of a total derivative 

4 E.B = = l -d»{e^ Xa A»d x A°) (6) 

where EM field tensor = d^A^—dvA^ and its dual is obtained using totally anti-symmetric fourth rank Levi-Civita 
tensor *F flv = ^€ l _ ll/a pF a P . Euler-Lagrange equations in this case have formal resemblance with axion electrodynamics 
of particle physics: Wilczek in 1987 [6] speculated that some of the axion field properties could possibly be realized in 
condensed matter system. This idea is nicely articulated for topological insulators in [2] arguing that observing elusive 
objects axions and monopoles in laboratory seems feasible. Thus TI and TMI have wider ramifications extending to 
field theory and particle physics. 

Attention is drawn to another significant idea in the unification quest for fundamental interactions that of duality 
[7]; it seems first work on duality in TI has been reported by Karch [8] in which properties of TI are discussed in 
the framework of SL(2, Z) symmetry. The present paper is aimed at understanding generalized duality in EM field 
equations such that the rotation angle is a local space and time dependent variable. In this formulation it is found 
that axion electrodynamics appears as a special case under suitable conditions and topological quantization is related 
with the multi-valuedness of what we term as duality gauge potential. This theory predicts a new effect that could 
be interpreted as dual to the magnetic field induced electric charge in the axion electrodynamics [6]. Local duality 
invariance gives new insights on the physics of TI as discussed in the following. 

That source-free Maxwell equations in vacuum have a duality symmetry E — > B and B — > — E has been known 
since long [7] . Postulating magnetic charge and current Maxwell equations in vacuum can be written in the covariant 
symmetric form 

8 V F^ = 4ttJ» (7) 



d v *F^ = 4ttJ£ (8) 
These equations are invariant under duality rotation with a real constant angle £ 

= F^cosC, +* F^sinC, (9) 



*F liV = -F^sinC +* F^cosC (10) 

provided the source current densities also undergo similar rotation. One of the important applications of duality 
symmetry is that one could explain the absence of magnetic charges [9]. Making a duality rotation with an angle 
such that the transformed = one gets the standard Maxwell equations: thus it is by convention that electron is 
assigned an electric charge — e and a vanishing magnetic charge. 

In a material medium duality rotation has to be applied carefully taking into account the constitutive relations (1) 
and (2) since the fields D and H would appear in the field equations [8]. Here we give a different approach introducing 
antisymmetric polarization tensor P^ and construct its dual *P> iV using the Levi-Civita tensor. The components P° l 
and P* J represent the electric P and magnetic M (P 4 - 7 = — e yfe Mfe) polarization vectors respectively. We can define 
the polarization current densities 

J£ = d v P v » (11) 



* J£ = d v *P^ (12) 
and modify the Maxwell equations in material medium in the presence of sources in the following form 

0„iW = 4nJ£ + 4nj£ = 4nJ>* (13) 



d v *F V * = 4ttJ^ + 4tt * J£ = 4tt * 



(14) 
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Obviously Eqs (13) and (14) are invariant under the duality rotation assuming that the polarization current densities 
also transform in the same way as the source current densities. Invoking the argument that one could choose the 
rotation angle such that the transformed * J M vanishes we obtain the Maxwell equations containing the macroscopic 
fields D and H in the traditional form. 

It is now straightforward to implement local duality invariance in Maxwell equations making £ to be an arbitrary 
function of space and time: a duality gauge 4- vector W v has to be introduced that transforms as W u — > W v + g~ 1 d I/ £, 
and couples to the EM field tensor in the field equations with the coupling constant g 

d v F v » = gW v *F V » + 4ttJ^ (15) 

d v = -gW v F v » + 4tt * J" (16) 

Generalized Maxwell field equations (15) and (16) constitute one of the important results of this paper. This set 
of equations in vacuum and source-free case was obtained earlier using variational principle for a vector Lagrangian 
[10], however there are subtle technical aspects in connection with such a formulation for local duality discussed in 
detail elsewhere (in preparation) . Here we work with the field equations and show that novel properties of topological 
insulators [2, 3] could be understood based on them, and put forward the proposition that phenomenological physical 
mechanism for new electric and magnetic properties of matter, e. g. TI and TMI is related with the duality invariance. 

Let us analyse the structure of Eqs (15) and (16): right hand sides comprise of three terms, namely duality induced 
term, polarization currents and true (free) charge currents. To understand them first we go back to Eqs. (13) and 
(14). In a usual way one can define new fields G^ v = F^ v — P^ v and its dual and rewrite these equations in the form 

d v G^ = 4nJ? (17) 

d v *G V » = 4ttJ£ (18) 

Alternatively since the polarization currents in a material originate due to bound charges i. e. electric dipole, magnetic 
dipole and multipoles, these could be clubbed together with the source currents as depicted in Eqs. (13) and (14). 
Thus we have two equivalent descriptions: modifying EM fields through constitutive relations as a material response as 
in Eqs. (17) and (18) and changing source currents in the presence of matter as in Eqs. (13) and (14). Extending this 
argument to Eqs. (15) and (16) the duality induced effect could be interpreted either in terms of altered constitutive 
relations or as a source of new current. 

To make the physical interpretation more transparent let us first specialize to the case when *Jp = 0, and assume 
simple relations (1) and (2) to write vector form of Eqs. (15) and (16) 

V.D = 5 W.B + 47rp e (19) 
V.B = -gW.E + 4irp m (20) 

V x E + — - 5 WxB| gW E = 4^J m (21) 

VxH-?+ 5 WxE| gW B = 47rJ e (22) 
at 

Conventionl assignment of discrete symmetry for electric charge density [9] implies that W and Wo are odd under 
both P and T transformations. We interpret as duality gauge field since it affects duality rotation, and seems to 
mediate spin interaction as reflected in the form of Eqs. (19)-(22). Further support to this interpretation could be 
given recalling the considerations on relativistic equations of motion for spin, see e. g. Section 11.11 in [9]: Thomas 
precession of spin can be understood using P and T odd 4-vector that couples to the EM field tensor F^ u in the 
form S V F^ U to give the equation of motion. Moreover note that discussing the physics of TMI [5] authors argue that 
static axion field could become dynamical in the antiferromagnetic phase due to spin wave and amplitude excitations 
underlining the significance of spin interaction in this phenomenon. 

Remarkably becomes axion- like field assuming W = — Va and Wq = — -|f , and setting magnetic charge current 
zero when (19) to (22) reduce to 

V.D = -gVa.B + 4np e (23) 
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V.B = (24) 

vxe + £- s |e = o < 25 » 

<9D da 
V x H - — - g(Va x E + — B) = 4ttJ e (26) 

To obtain above equations it is necessary that Va is normal to E and parallel to B: this is an important restriction that 
does not follow from the derivation of axion electrodynamics using the topological action (3) in [2] . Further Faraday's 
law (25) also gets modified when a is time-dependent as compared to axion electrodynamics in which Faraday's law 
is unaltered. What about the key ingredient: topological quantization? Though there is no topological action in 
the present formulation topological quantization arises if duality gauge field is not single- valued. It is similar to the 
topology of Aharonov-Bohm effect [11]: in the standard illustration a multiply connected space follows for a — <j> in 
cylindrical coordinates (r, <f>, z) in which periodic boundary <fi — > <j> + 2n changes the field a. Once this topological 
quantization is used the arguments given by Qi, Hughes and Zhang in [2] could be immediately adopted to arrive at 
the new properties of TI i. e. TME and image monopole, since the field equations (23) to (26) are identical to their 
equation (104) for static a. Interestingly the gauge coupling constant may also be obtained using multi-valuedness 
of a and transforming <?Va.B in (23) to g4>&o where $o = ff is magnetic flux quantum: changing a from to 2n 
induces an electric charge, if it is assumed to be 47re then we get g — 2a/ n. 

Preceding discussion shows that duality invariant theory explains all the novel features of topological insulators for 
static a, and predicts a modification in Eq. (25) for dynamical a. Expected consequence of modification in Eq. (25) 
would be that the dispersion relation of so called axionic polariton proposed in [5] will change. Next we discuss an 
entirely new effect that emerges in the duality invariant theory. 

Let us consider dual of axion-like electrodynamics: assume Va normal to B and parallel to E then Eqs. (19) to 
(22) reduce to the following equations 

V.D = 47rpe (27) 
V.B = .aVa.E (28) 
da 

V x E + — + g(Va x B + — E) = (29) 
<9D Bo 

VxH-^-.g^B = 47rJ e (30) 

Eq. (28) shows that electric field induces a magnetic charge, and corresponding magnetic current in Eq. (29). 
In contrast to axion electrodynamics its dual is thus seen to be characterized by a purely duality (topological) 
electrodynamics as free magnetic source current and material parameters e and /U are not present in Eqs. (28) and 
(29). Note that axion field effect in TI is relatively very small as compared to free current J£ and the parameters 
e and \i and though one is looking for surface effects their observation becomes very difficult. Therefore search for 
possible materials and viable experiments to detect the dual effect embodied in Eqs. (27) to (30) would be of great 
value. 

An intriguing aspect regarding induced magnetic charge is worth mentioning. To calculate it in Eq. (28) we perform 
the similar steps as we did for calculating induced electric charge in Eq. (23). Assuming electric flux equal to 47re 
induced magnetic charge is obtained to be 167rea. Recalling that monopole charge according to Dirac quantization 
condition is e/2a, the ratio of induced magnetic charge to monopole charge turns out to be extremely small: 327ra 2 . 

In the light of symmetry and universality for TI electrodynamics emphasized in [2] it becomes natural and attractive 
to explore generalized duality invariance. However it is imperative to seek its foundation or link with some variant of 
band theory that takes into consideration explicit calculation of various interactions in the material. Dynamical axion 
in MTI may have origin in the coupling of electrons with the antifcrromagnctic order parameter [6] , therefore it would 
be interesting to examine the role of charge density wave order parameter afresh in conjunction with ferromagnetic 
and antiferromagnetic order parameters in new materials if the dual to axion field predicted in the present work could 
be realized. In fact, discussing SL(2, Z) symmetry in the case of dynamical axion field Karch [8] has rightly noted that 
duality demands e and \i also to become dynamical fields - the dilaton. In our duality invariant theory, the complete 



■5 



set of Eqs. (15) and (16) shows that the material parameters e and /x need not be dynamical variables rather the 
dynamical axion field is a manifestation of the duality gauge field. 

Before we conclude, the significance of topological phases in optics, in particular, Pancharatnam phase for unraveling 
the unique properties of TI, for example, TME is pointed out. Though optical techniques have been extensively used 
to study interfaces and surface properties of insulators a recent experiment [12] seems important in which TI is probed 
using ultrafast laser pulses in the nonlinear regime. Duality invariant rendition of TI electrodynamics presented here 
indicates the role of spin angular momentum and topological phase of polarized light besides the topological property 
of the insulating material. We suggest that the kind of experiment performed by Hsieh et al [12] could be enlarged 
in scope to search for any such topological phase effect in TI. One possibility would be to study optical rotation of 
polarized light incident on vacuum - TI - ordinary insulator and observe reflected and transmitted light for optical 
rotation [2, 8], if any, and analyze the observations for Pancharatnman phase. 

In conclusion a generalized duality invariant electrodynamics is developed applicable to topological insulators and 
it is shown that under appropriate conditions this theory reduces to the axion field-like theory widely discussed in the 
literature; a new effect dual to axion electrodynamics is also predicted. 
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